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Abstract 

Within a dynamical transport approach we investigate charged current interactions 
in neutrino-nucleus reactions for neutrino energies of 0.3 - 1.5 GeV with particular 
emphasis on resonant pion production channels via the A33(1232) resonance. The 
final-state-interactions of the resonance as well as of the emitted pions are calculated 
explicitly for 12 C and 56 Fe nuclei and show a dominance of pion suppression at mo- 
menta p n > 0.2 GeV/c. A comparison to integrated ir + spectra for f^+ 12 C reactions 
with the available (preliminary) data demonstrates a reasonable agreement. 
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1 Introduction 

The interest in neutrino induced reactions has been rapidly increasing during the 
last decade while current and future long base-line (LBL) experiments, K2K, Mini- 
Boone, MINOS, CERN-GS, J2K etc. aim at determining the electro-weak response 
of nucleons and nuclei. In general the neutrino- nucleon/nucleus interaction consists 
of the quasi-elastic scattering, resonance production and deep inelastic scattering, 
the latter becoming dominant only for neutrino energies above about 1 GeV increas- 
ing further linearly with energy. Whereas at low neutrino energies E v < 0.3 GeV the 
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quasi-elastic channel clearly dominates, one observes an increasing contribution from 
resonance excitations at higher energies up to about 2 GeV [1] . Here charged current 
(CC) reactions offer the unique possibility to measure the /i~ in coincidence with e.g. 
the charged 7r ± or neutral it in order to extract information on the resonance tran- 
sition form factors since the 4-momentum of the lepton fixes the transferred energy 
v and 4-momentum squared q 2 . 



A variety of models (and event generators) treat neutrino nucleus reactions [2,3,4,5,6,7,8] 
and incorporate selected aspects of scattering on nuclei like Pauli-blocking, nuclear 
Fermi motion, nucleon binding energies, pion rescattering and charge exchange reac- 
tions. However, most of the models lack an explicit propagation of the resonances and 
their in-medium decay to pions and nucleons as well as more realistic coordinate (and 
momentum) distributions of the target nucleons. In order to study the importance of 
these effects we employ in this work a transport approach that has been tested in a 
wide range of n, p reactions on nuclei as well as relativistic nucleus-nucleus collisions. 
For reviews we refer the reader to Refs. [9,10]. This transport model - or related 
versions - have been also used and tested in the description of the photo-nuclear 
and electron induced reactions with the momentum and energy transfer compara- 
ble to those considered here [11,12] as well as at much higher energies [13]. This is 
particularly relevant, since the nuclear response to electromagnetic probes is quite 
analogous to the weak probes addressed in this study. 



For the description of the elementary process on nucleons we use the model of Rein 
and Sehgal [14] - with slight modifications - that has been often employed for an 
estimate of neutrino cross sections for resonance transitions. It is based on quark 
oscillator wave functions to evaluate the form factors as developed by Feynman 
et al. [15]. Alternatively, one may also parametrize the neutrino-nucleon-resonance 
vertex by introducing 'proper' vector and axial-vector form factors. The differential 
cross sections are well known [16] and usually given in terms of helicity amplitudes. 
Our present investigation is focused on final-state-interactions (FSI) of the excited 
A(1232) resonance and the emitted pion. Total ir + cross sections from 12 C targets 
will be presented in comparison to preliminary data as a function of the neutrino 
energy E v from 0.3 to 1.5 GeV. We also provide energy differential 71^,71° spectra 
for 12 C and 56 Fe at E v — 1 GeV and perform a differential analysis of final-state- 
interactions (FSI) on both targets that are used in the MiniBoone and MINOS 
detectors as target material. Note that alternative forms of the neutrino-nucleon- 
resonance vertex will modify the double differential cross sections but have no severe 
impact on the strength of the FSI. 
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2 Neutrino-nucleon reactions 



We consider the interactions of a neutrino with a nucleon that lead to a muon 
and pion in the final state, i.e. the charged current (CC) interactions 

Vp, + V — > /J>~ + 7T + + p (1) 

+ n — > fx~ + 7r + + n (2) 
z/ M + n — > fi~ + 7T° + p. (3) 

Neutral current (NC) interactions are 

Vp+P^Vn + iP+p (4) 

Vft+P ^ + + n (5) 

v» + n -> z/ M + n~ + p (6) 

+ n -> z/ M + 7T° + n. (7) 

The NC reactions with the 7r° production, eqs. (4) and (7), represent a dangerous 
background in the long-baseline neutrino oscillation studies, since they could be 
confused with the quasi-elastic electron production caused by neutrino-oscillations. 
We will consider the NC pion production separately in future work. 

We denote the energy transfer by v — E v — Ei, where E\ is the energy of the final 
lepton, the 4-momentum transfer squared by q 2 = v 2 — Q 2 = M 2 — m 2 n — 2m n v with 
M,m n ,Q denoting the resonance mass, nucleon mass and 3-momentum transfer, 
respectively. We then obtain for the invariant energy (with the nucleon at rest) 
W = v + m n or alternatively W 2 = M 2 + Q 2 . The double differential cross section 
for neutrino A-resonance excitation then reads 



da (G F cos(0 c )) 2 q 2 



dq 2 du 47r 2 Q 2 

with 



k [u <jl + v or + 2uva s ] (8) 



E„ + Ei + Q E u + Ei - Q M 2 -m 2 n 

u = , v = , k = . 9 

3E V 3E U 2m n y J 



The a terms read 

^,*(g 2 ,W0 = — O^Tll^lA,^)! 2 2WA(W,M), (10) 
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Wm O 2 

a s (q\W) = -^^Y.\( N >3^1\F \A, Jz )\ 2 2WA(W,M), (11) 



kM q 1 



3z 



where A(W, M) denotes the A spectral function - to be specified below - which 
replaces the conventional term 2ttW5(W 2 — M 2 ). The N — A transition amplitudes 
F T , F in (10, 11) are well known in principle [14], however, involve rather uncertain 
vector and axial-vector form factors which we adopt in the form 



G V ' A (q 2 ) 



,2 / n 2 



-2 



f-^kV-wJ (12) 



with M v ,a specifying characteristic mass (or inverse length) scales. We use the pa- 
rameters [17] My = 0.843 GeV and M A = 1.032 GeV in the following calculations 
(default parameter set A). 

The A spectral function is taken in the form (with Ma = 1.232 GeV) 

WT(W) 

A W M ^ = ( w>-Mi? + w^(wy (13) 

where the width Y explicitly depends on W . In vacuum a good parametrization is 
given by [18] 



with 

q 2 (W) = (W 2 - (m n + m^){W - (m n - m 7r ) 2 )/(4^ 2 ). (15) 

Following [18] the parameters are T = 110 MeV, (3 = 0.3 GeV. Furthermore, 
q r = q(M A ) while denotes the pion mass. We have checked that alternative 
parametrizations only moderately change the shape of the cross section (8) in line 
with Ref. [19]. 

The results for da/dq 2 dv for the channel +p — > p<r + A ++ are displayed in Fig. 1 
for a neutrino energy E v — \ GeV and show a characteristic bump in the transferred 
energy v due to the spectral function (13) which softens and kinematically shifts 
with increasing negative q 2 . A similar functional dependence on v and q 2 is known 
from e~ + A reactions [12] . 

Since presently no double differential data are available for a detailed comparison 
and fixing of the form factors (12), we use cross sections that are integrated over v 
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Fig. 1. The double differential cross section da/dvdq 2 for the reaction v + p — > \i~ + A ++ 
as a function of the transferred energy v and the 4-momentum transfer squared q 2 for a 
neutrino energy E u = 1 GeV. 



and g 2 . Note that the data in Fig. 2 [20] correspond to reactions on if 2 and D2 and 
differ significantly from each other. Our standard 'historical' parameter set A [17] 
falls slightly low for both reactions, while calculations with a modified parameter 
set B {My = 1 GeV, M A = 1.2 GeV) perform better (though not perfectly). Our 
model gives a ratio of the 7r + cross section from the proton relative to the neutron 
of 3 which is determined by a relative Clebsch-Gordon coefficient of v3 in the decay 
amplitude. The original Rein-Sehgal model [14] has an additional factor \/3 in the 
production amplitude for the A ++ relative to the A + , which gives a total ratio of 1/9 
for the cross sections + n — ► /i - + n + + n relative to + p — ► ^ + n + + p. Since 
the data in Fig. 2 are not compatible with a ratio 1/9 but favor a ratio 1/3 we have 
'redefined' the production matrix element to be identical for A ++ and A + . Note that 
also higher baryon resonances may contribute to 7r + production - especially for the 
Vn + n — > + 7r + + n channel - which are discarded in this study. These resonances 
have independent transition form factors that presently are not controlled by data 
and require further parameters. The explicit dependence of the cross section on E v 
reflects both the form factors employed as well as the detailed form of the spectral 
function (13). 

In view of the large uncertainties in the data we do not aim at fitting the form 
factors here and consider instead reactions on nuclei by employing a coupled-channel 
transport model. 
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Fig. 2. The calculated cross sections for the reactions v + p — > \i + 7T + + p (upper part) 
and v + n — * fj,~ + 7r + + n (lower part) as a function of the neutrino energy E v for the 
standard parameter set A (M v = 0.843 GeV, M A = 1.032 GeV) and the modified set B 
{My = 1 GeV, Ma = 1.2 GeV). The data are taken from Refs. [20], respectively. 

3 Neutrino nucleus reactions 



Relativistic transport models have been originally developed for nucleus-nucleus col- 
lisions starting from the Fermi-energy domain up to ultra-relativistic energies (y/s ~ 
200 GeV) in order to extract physical information from nucleus-nucleus collisions 
relative to scaled proton-nucleus reactions [9,10,21,22,23,24]. These models incorpo- 
rate a variety of in-medium effects such as mean-field potentials - i.e. the real part of 
the retarded selfenergy - as well as finite life time effects of hadrons due to multiple 
interactions (decays) in the medium. The latter are encoded in dynamical spectral 
functions for the hadrons of interest that are fully specified in terms of the real and 



6 



imaginary parts of the hadron selfenergy. We recall that an off-shell transport version 
- suitable for the propagation of particles with finite spectral width - is available, 
too [25]. While the initial studies had been devoted to A + A collisions, these models 
have been used in the description of pion-nucleus [26] and photo-nuclear reactions 
[11,12,27] in more recent years as well. The advantage of such transport models is that 
they are suited for a wide variety of reactions and thus can relate e.g. photo-nuclear 
reactions to pion/proton induced reactions as well as neutrino-nucleus reactions (as 
addressed here). 

In this work we employ the CBUU 1 approach [28] which is based on a set of transport 
equations for the phase-space distributions fh(x,p) of hadron h that are coupled via 
collision terms describing the mutual interaction (or decay and formation) rates in 
phase space. The mean-fields for the baryons are the same as in Ref. [24] while the 
description of hadronic collisions is identical to Ref. [10] . 

For the neutrino-nucleus reactions we employ a perturbative treatment as in the 
case of photo-nuclear reactions [13] which allows us to obtain good statistics even 
for weakly interacting probes. In this respect the total neutrino-nucleus reaction is 
divided into two parts: i) the initial neutrino interaction on a single nucleon of the 
target (including its local Fermi motion and mean-field potential) and ii) the prop- 
agation and decay of the excited resonance (as well as the decay products) in the 
residual nucleus. Each excited resonance in step i) is assigned a weight Wi propor- 
tional to the double differential cross section da/dudq 2 where v and q 2 are selected by 
Monte Carlo in the kinematical allowed regime according to the differential distribu- 
tions displayed e.g. in Fig. 1 for E v = 1 GeV. A large number of simulations (~ 10 4 ) 
are performed such that the sum over the weights Wi - normalized by the number 
of events - is equal to the v and q 2 integrated cross section per nucleon. Note that 
for fixed v and q 2 and each individual 4-momentum q^ of the nucleon the invariant 
energy as well as the 3-momentum of the excited resonance are fully determined by 
the energy and momentum conservation. 

We then follow the motion of the perturbative hadrons (here A's, decay pions and 
nucleons) within the full background of residual nucleons by propagating the A's 
with 2/3 of the nucleon potential - neglecting pion potentials in this study - and 
compute their collisional history with the residual nucleons of the target. A reduced 
A potential - relative to the nucleon potential - is adopted in order to improve 
the description of pion spectra from p+A reactions and to incorporate the findings 
from Peters et al. [29] for coherent photo-production of pions on nuclei. In order to 
trace back the FSI of the hadrons the A-production vertices are recorded with their 
individual weight (including their 4-momenta) as well as the 4-momenta of the final 
pions. 

We point out that an alternative transport description of neutrino-nucleus reactions 
on the basis of the GiBUU transport model has been initiated by T. Leitner et al. 

1 Coupled-Channel Boltzmann-Uehling-Uhlenbeck 
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in Refs. [30]. 



Here we present the calculated cross sections for ir + meson production in CC reactions 
on 12 C targets. The results are displayed in Fig. 3 for the two parameter sets discussed 
in the context of the elementary cross sections shown in Fig. 2. We find that the 
parameter set A underestimates the preliminary data on 12 C from Ref. [31] as in 
case of the proton or neutron target in Fig. 2 while the parameter set B performs 
better (as expected from Fig. 2). Note, however, that any agreement /disagreement 
is primarily due to the 'uncertain' axial- vector form factor G A (q 2 ) such that no final 
conclusions can be drawn from the comparison in Fig. 3. 




0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 



E (GeV) 

Fig. 3. The total cross section a n + for the reaction v + 12 C — > [i~ + ir + + X as a function of 
the neutrino energy E v for the two parameter sets discussed in the text. The preliminary 
data have been taken from Ref. [31]. 

As an example the differential pion (tt + , 7r°) spectra in the pion kinetic energy T n - 
divided by the pion laboratory momentum Pj a & squared - are shown in Fig. 4 for a 
12 C and 56 Fe target at the neutrino energy E u = 1 GeV. The differential spectra 
da^/dT^/ P 2 ab for 12 C and 56 Fe show a shoulder at low T n and are roughly exponential 
for high pion energies. The shape of the spectra is similar for both targets but they 
do not simply scale by a factor of 36/8 = 4.5 as expected from a scaling of the 
elementary cross sections. (This factor comes about as follows: (26+30/3)/ (6+6/3) 
when counting the protons and neutrons with their excitation to A ++ and A + , 
respectively, and the relative decay to tt + . For tt° the decay fraction of the A + 
is 2/3, however, no direct 7r°'s can be produced in CC reactions on protons. The 
corresponding scaling factor thus is 30/6= 5.) In case of the Fe target we also display 
the 7T~ spectra (upper part) that are down by an order of magnitude relative to the 
7r + spectra and stronger localized at low pion energies. Note that no direct negative 
pions can be produced in the CC neutrino reactions on nucleons and therefore the 
final 7r~ (or tt~/tt + ratio) provide a sensitive probe for the strength of the FSI. 
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Fig. 4. The differential cross sections da n + /dT w / Pf ab (upper part) and da^o / dT n / Pj* ab 
(lower part) for the reactions v + 12 C(( 56 i ? e) — * \i~ + tt + X as a function of the pion 
kinetic energy T n in the laboratory for a neutrino energy E u = 1 GeV. The tt~ differential 
spectrum (in case of the 5e Fe target) is shown by the dashed line (in the upper part) and 
is entirely due to final state interactions in the target nucleus. 

Since the differential pion spectra suffer from the same uncertainties as the integrated 
spectra we turn to ratios of cross sections that reflect dominantly the influence of 
the FSI. For this purpose we consider the ratios 

_ (d 3 a x /d 3 p) fin 

RM ~ W^Whn (16) 



where the initial differential spectra in the denominator of Eq. (16) are determined 
from the decay of the primarily produced A resonances without FSI. Fig. 5 shows 
the ratio Rq of the final differential ir° momentum spectrum to the neutral pion 
spectrum without FSI for + 56 Fe at E v = 1 GeV as a function of the longitudinal 
momentum P z and transverse momentum Pt (with respect to the direction of the 
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incoming neutrino). This ratio is roughly constant (P ~ 0.4) for yP^ + P 2 > 0.2 
GeV/c and exceeds unity for lower pion momenta in the laboratory due to strong 
FSI of the A's and decaying pions with nucleons. We recall that NA interactions in 
the medium involve - apart from energy and momentum exchange - charge transfer 
as well as AiV — > NN absorption reactions. The latter cross sections are obtained 
from refined detailed balance relations that incorporate the finite (mass dependent) 
width of the A- resonance [28]. The strong reduction of the pion spectrum at higher 
pion momenta is essentially due to 7r + iV — > A resonant scattering. 




Fig. 5. The ratio Rq of the final differential n momentum spectrum to the it spectrum 
without FSI for Vn + 56 Fe at E u = 1 GeV as a function of the longitudinal momentum P z 
and transverse momentum (with respect to the direction of the incoming neutrino). 

Apart from the overview in Fig. 5 it is, furthermore, of interest to quantify the 
strength of the FSI e.g. as a function of the pion kinetic energy in the laboratory 
T n . All ratios displayed in Fig. 6 for tt + (upper part) and tt° (lower part) show a 
characteristic dip at low kinetic energy T n followed by a bump (with a ratio above 
unity due to pion rescattering) and an approximate plateau for T n > 0.15 GeV which 
is lower for the heavier Pe-target than for the light C-nucleus. These differential FSI 
effects can also be directly assessed by taking experimental differential ratios of 
spectra for heavy and light targets. 

Note that charged current reactions with pion production and subsequent absorption 
of the pion might be misinterpreted experimentally as the quasi-elastic CC reactions 
and lead to an overestimation of the CC quasi-elastic cross sections. In case of the 
12 C target we find an integrated n + (ir°) absorption of 12% (14%). When impos- 
ing a detection threshold for pions of T n = 100 MeV in the laboratory the relative 
absorption amounts to 33% (35%) for 7r + (tt°). For the 56 Fe target the respective 
numbers are 35% (42%) for the integrated spectra and 58% (64%) for an energy cut 
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Fig. 6. The differential ratios R + (T W ) (upper part) and Rq(T w ) (lower part) for the reac- 
tions Un + 12 C( 56 -Fe) — > \i~ + 7T + (7r°) + X as a function of the pion kinetic energy T T in 
the laboratory for a neutrino energy E v = 1 GeV. The error bars denote the statistical 
uncertainty in determining the ratios numerically due to the low cross sections at large 
pion energies (cf. Fig. 4). 

4 Summary 

In this study we have reported first dynamical transport calculations for neutrino- 
nucleus reactions in the energy range from 0.3 - 1.5 GeV. The advantage of our 
transport description is that it is suited also for a wide variety of reactions and 
thus can relate neutrino- nucleus reactions to photo-nuclear reactions or pion/proton 
induced reactions on nuclei. Furthermore, standard in-medium effects such as mean- 
field potentials or in-medium broadening of hadronic resonances can be exploited 
within the same framework. Effects of Fermi motion and Pauli-blocking of nucleons 
are incorporated by default. 

The method employed here is quite general and can be applied to a wide variety 
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of nuclei and reaction channels that incorporate all possible final state interactions. 
Furthermore, the transport approach can be used for an event-by-event analysis and 
correlate the momenta of different particles in each event, e.g. emitted nucleons with 
pions etc.. It also allows for a clean implementation of experimental acceptance cuts 
(i.e. detector thresholds and/or detector geometry). Hence, it should be applicable 
to the description of the neutrino-pion production on nuclei, in particular to the 
description of the final-state-interactions (FSI) of the produced hadrons. 

In our first investigation we have focused on resonant pion production channels 
via the A33(1232) resonance in charged current (CC) reactions. The final-state- 
interactions of the resonance as well as the emitted pions have been calculated ex- 
plicitly for 12 C and 56 Fe nuclei and show a dominance of pion suppression for pion 
momenta above 0.2 GeV/c due to A absorption on nucleons (iVA — > NN) as well 
as strong pion rescattering in the kinematic regime of the A resonance excitation. 
Neutrino reactions on 16 are similar to 12 C targets and will be reported elsewhere. 
A comparison to integrated 7r + spectra for + 12 C reactions demonstrates a rea- 
sonable agreement in view of the uncertainties in the vector and axial-vector form 
factors for the elementary reactions on nucleons. 

We stress that precise double differential data (with respect to v and q 2 ) on free 
nucleons are mandatory before addressing the FSI in a differential manner as pro- 
posed in this study. Once this task is solved the question of in-medium effects on 
resonance transitions in neutrino-nucleus reactions can be addressed and transition 
amplitudes to higher baryonic resonances can be investigated. The coupled-channel 
transport calculations used here already include the various in-medium phenomena 
and provide a firm basis for the analysis of future differential data. 

The authors like to acknowledge discussions with T. Leitner who published a related 
study in the GiBUU model [32] after submission of our manuscript and pointed out 
an error in our earlier calculations. 
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